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SUMMARY 


A method is described herein for the prediction of rolling moments on a wing 
penetrating a vortex velocity field generated by a large aircraft. Rolling moments 
are determined from lifting pressure coefficients computed with an inviscid-flow 
linear panel method. Two empirical corrections are included to account for the 
lifting efficiency of an airfoil section and the local stall on the wing. Predicted 
rolling moments are compared with those from two wind-tunnel experiments. Results 
indicate that experimental rolling moments, for which the Reynolds number of the 
following wing is low, should be interpreted with caution. 


INTRODUCTION 

One of the greatest hazards to a small aircraft in an airport environment is a 
strong rolling moment induced by an axial penetration of a wake vortex produced by 
a large aircraft. Determination of the induced rolling moment is, therefore, a 
vital part of the study of the wake vortex hazard. The effectiveness of a given 
vortex attenuator can be studied by measuring the vortex velocity distributions in 
the wake of an aircraft with and without the attenuator and calculating the relative 
rolling moments produced on a following aircraft. 

The validation of a predictive method requires measured wake velocity distri- 
butions as well as induced rolling moments on a wing in the same wake. However, 
only a limited amount of data is available which includes both measurements. 
References 1 to 3 present this type of data as well as some comparisons with 
applicable theories. Predicted rolling-moment results from two-dimensional strip 
theory, strip theory with an empirical lift-curve slope correction, and vortex- 
lattice theory are presented in reference 1. In reference 2 rolling moments are 
predicted with simple strip theory in which the local lift-curve slope is assumed 
to be the same as the overall measured lift-curve slope and the lift is limited to 
a maximum value. Inviscid-flow rolling-moment results from lifting-line theory are 
presented in reference 3. The present approach (with a linear panel method) involves 
corrections based on measured airfoil data for both local lift-curve slope and local 
stall (maximum lift) . 

The purpose of this paper is to present a method for computing the rolling 
moment for an axial penetration of a wing into a vortex whose velocity distribution 
is known. Since the method accounts for variations with Reynolds number, it can 
be useful in the interpretation of wind-tunnel data. Development of the method 
was stimulated by work in the Langley Vortex Research Facility (refs. 4 and 5) and 
by the introduction of the laser-velocimeter technique (ref. 6) for measurement of 
velocity distributions. Results in the present paper are compared with experimental 
data from references 2 and 3. 



SYMBOLS 


a area of finite-element panel, in 

B span of vortex-generating wing, in. 

b span of following wing, in. 

rolling-moment coefficient (eq. (1)) 

C lift coefficient of vortex-generating wing 

c chord of following wing, in. 

c^ section lift coefficient 

c. section lift-curve slope, per deg 

F section lift factor (eq. (7)) 

R c Reynolds number based on chord of following wing 

2 

S following-wing planform area, in 

V/U^ ratio of vortex velocity to free-stream velocity 

x rectangular coordinate in streamwise direction, in. 

y rectangular coordinate in spanwise direction, in. 

a angle of attack, deg 

a effective stall angle (eq. (8)), deg 

ACp lifting pressure coefficient on a finite-element panel 

Subscripts : 

exp experimental 

i summation index on wing panels 

L left group of finite-element panels 

max maximum 

o centerline of following wing 

R right group of finite-element panels 

A prime on a symbol indicates the value is calculated from inviscid-f low linear 
theory for a flat-plate wing? a bar over a symbol indicates a coordinate at the 
centroid of area of the finite-element panel. 
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METHOD 


A vortex velocity field imposes a torque which causes an aircraft encountering 
the vortex to roll. For an axial penetration, the initial rolling moment is the 
maximum value. Throughout the present paper, the term "rolling moment" means 
"initial rolling moment." 


The rolling moment on a wing in an axial encounter with a wake vortex can be 

approximated fairly well by inviscid-f low steady-state theory (ref. 7) , but only for 

very large Reynolds numbers . Inviscid-f low calculations are not sufficient for 

predicting the rolling moment in a wind tunnel where the Reynolds number of the 

following wing is very low. The steady-state method described herein includes 

empirical corrections to account for the effects of Reynolds number and of airfoil 

camber and thickness on the lift an airfoil section can attain. A computer program 

by Carmichael et al. (ref. 8) for inviscid linear-theory flow over wing-body 

combinations is used to calculate the pressure distribution on a flat-plate wing 

which is modeled by a number of finite-element panels. An input modification allows 

for the nonuniformity (vortex velocity distribution) in the incoming airstream and 

for local stall (c, ) on the wing. A section lift factor is used to account for 

i , max 


the experimental lifting efficiency (c. ) of the airfoil. The resulting wing pres- 

l a 


sure distribution is used in the computation of the rolling-moment coefficient. 


The effect of the following wing on the vortex velocity distribution is not 
considered in the present work. The fuselage is not modeled and Mach number effects 
are not taken into account. Also, post-stall lift loss is not considered in the 
local stall correction. 


Wake Vortex Simulation 

The effect of a nonsymmetrical vortex velocity distribution (e.g., fig. 1) on a 
following wing is accounted for by imposing an artificial twist distribution on the 
wing. The local wing incidence is the sum of the angle of attack a, the geometric 
twist angle, and the artificial twist angle (arctan (V/Uoo) ) • To account for stall, 
the local wing incidence is limited by the angle corresponding to c ^ /Inax for the 
test Reynolds number based on the following-wing chord. 

The linear- theory aerodynamics program (ref. 8) models a wing with a number of 
constant-pressure panels. It requires the geometry of a semispan wing as input 
and automatically reflects across a plane of symmetry to account for the other half 
of the wing. (See example in fig. 2(a).) Therefore, any twist distribution 
specified by program inputs appears symmetrically in the reflected geometry. The 
appropriate nonsymmetrical twist distribution cannot be imposed on the wing in the 
usual way. Therefore, the input wing is taken as the entire wing so that a 
nonsymmetrical twist distribution can be applied as needed. The aerodynamic inter- 
ference from the reflected wing must be reduced by moving the input wing away from 
the program symmetry plane (fig. 2(b)). At a separation distance of 10 wing spans 
(fig. 2(c)) the interference between the input and reflected wings results in less 
than 0.1 percent error in the lifting pressure coefficients. A separation distance 
of 10 or more wing spans is used in practice to ensure negligible interference 
between the wing of interest and its reflected image. A large spearation distance 
is achieved by simply adjusting inputs to the program. 
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Rolling-Moment Coefficient 


The rolling-moment coefficient is computed from the lifting pressure coefficient 
ACp and the spanwise rectangular coordinate of the centroid of area y for each of 
the 100 finite-element panels that make up the input wing. Five panels are used in 
the chordwise direction and numbered as shown in figure 3(a). Since the input wing 
is taken as the entire wing, the panels are renumbered into a left group and a right 
group, as shown in figure 3(b). Thus, 

C \ = c l ,L + C l ,R d) 


The contribution from the left group is 
-i 50 


l,L - ib^ % l L,i a i,i' 5 L,i 


Similarly, 


-i 50 

■o = — y ' (Ac ) .a (y 
,R sb ^ p R,i R, l ^R, i 


“ y. 


(2) 


(3) 


The lifting pressure coefficients are determined from the section lift factor F 
(which is explained in the next section) and lifting pressure coefficients calcu- 
lated from inviscid-f low linear theory on a flat-plate wing: 


(ACp) 


L,i 


F (AC ) . 

P L, 1 


(4) 


and 


(Ac )_ . = f(Ac' )„ . 
p R, l p R, l 


(5) 


Since the constant F is a factor in all terms of equation (1) , 


C 1 " FC 1 


(6) 


Empirical Corrections 

Two empirical corrections are included in the rolling-moment coefficient; one 
corrects the airfoil section lift-curve slope and the other allows for local stall 
on the wing. The two corresponding parameters are the section lift factor F and 
the effective stall angle a es . These parameters can be determined from curves for 
airfoil experimental c as a function of a from sources such as references 9 
to 11- 

Section lift factor.- The section lift factor F is the ratio of (c, ) to 

1 a exp 

c* , where c' = 0.094 is calculated at the center of a rectangular wing of aspect 
ratio 1000 with no twist. Therefore, 
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F = (c ) /0.094 

La exp 


( 7 ) 


and accounts for viscosity and for airfoil camber and thickness. Consider an 

NACA 0012 airfoil at a low Reynolds number. The plot of experimental c^ as a 

function of a is shown in figure 4 for this airfoil at the lowest test Reynolds 

number (R c = 170 000) for which NACA 0012 data are given in reference 9. Note 

that the curve is linear with (c. ) = 0.100 below a = 5°. Therefore, 

La exp 

F = 0.100/0.094 =1.06 in this example. This factor corrects the calculated c^ 
so that it agrees with experiment at low angles of attack and follows the line of 
long dashes in figure 4. 


Effective stall angle.- The effective stall angle is 


a 

es 


(c ) /(c 

l,max exp 1 


) 

a exp 


(8) 


For the rolling-moment calculation, a es is used to limit the calculated c^ to 

the experimental maximum. For the curve in figure 4, (c, ) = 0.83 as shown 

1 ,max exp 

by the line of short dashes. Thus, from equation (8), a QS = 0.83/0.100 = 8.3°. 

(See intersection of dashed lines.) Since the airfoil is symmetric, a es = -8.3° 
for negative angles of attack. For a cambered airfoil, a different a es would have 
to be determined for negative angles of attack. To account for stall on the wing 
penetrating the vortex, the local wing incidence (geometric plus vortex-induced 
incidence) is limited to the effective stall angle. Consider a penetrating 
(following) wing with no geometric incidence (a = 0 and no twist) , so that the 
local wing incidence is simply arctan (V/Uco) . For the vortex velocity distribution 
from figure 1 and a symmetric airfoil with a es = 10° , figure 5 shows how the 
tangent of a es limits the vortex velocity distribution. 


RESULTS AND DISCUSSION 

Two cases are treated based on availability of relevant experimental data. 

Both are for a central penetration of a vortex by a wing with a - 0, no geometric 
twist, and a symmetric airfoil section. In each case, two Reynolds numbers are 
considered? one is higher than the experimental value and the other is near the 
experimental value. The effects of F and a es are studied by computing rolling- 
moment coefficients with and without these empirical corrections. Predicted results 
are presented with no corrections, with one correction, and with both corrections, 
and results are compared with experimental data. As noted earlier, verification of 
the present method requires experimental data for which rolling-moment coefficients 
and corresponding vortex velocity distributions are known. The effect of the 
following wing on the vortex velocity distribution is not considered. Airfoil- 
section data of the following wing must also be known or calculable. Results are 
computed for one case from reference 2 and for one case from reference 3. 


Case A 

Measured velocity distributions are given in reference 2 for trailing vortices 
at x/B = 2.5 and x/B = 5.0 behind a 35° sweptback wing with B = 42.00 in. for 
lift coefficients of = 0.36 and C L = 0.74. The distribution with the highest 
maximum velocity, namely, that for x/B = 2.5 and = 0.74, is used for case A. 
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(See fig. 1.) The rectangular following wing with an aspect ratio of 7.5 has an 
NACA 642-015 airfoil section with roughness to fix transition, b = 10.00 in., 
c = 1.33 in., and R c * 110 000 in the experiment. 

For R n = 6 000 000 and the NACA 64o-015 airfoil, ( c . ) = 0.110 and 

° exp 

(c, ) =1.10 from the standard- roughness lift curve in reference 10. From 

l, max exp 

equations (7) and (8), F = 0.110/0.094 = 1.17 and a es = 1.10/0.110 = 10°. 

(a es = -10° for negative a.) As shown in figure 6, F increases the uncorrected 
rolling-moment coefficient from 0.070 to 0.082 and a es decreases it from 0.082 to 
0.080 for R c = 6 000 000. (Compare figs. 1 and 5 to see the cause for this small 
decrease . ) 

Since no experimental lift curve is known to be available for the NACA 642 _ 015 

airfoil with roughness for a Reynolds number near R c = 110 000, corresponding 

estimates of (c, ) and (c, ) are now made. The ratio of the maximum 

la exp If max exp 

lift coefficients at two Reynolds numbers is assumed to be the same for rough or 

smooth symmetric airfoils with the same thickness-to-chord ratio. Lift curves are 

available in reference 10 for the NACA 642-015 airfoil with and without roughness 

for R c = 6 000 000 and without roughness for R c = 3 000 000. As shown in this 

reference, (c, ) =1.33 for R_ = 3 000 000 and (c, ) = 1.45 for 

1 ,max exp c l, max exp 

R n = 6 000 000 without roughness. Therefore, (c, ) sr 1.10(1.33/1.45) = 1.01 

c . 1 ,max exp 

for a rough NACA 642“015 airfoil at R c = 3 000 000. For another smooth 15-percent- 

thick symmetric airfoil., namely, an NACA 0015 (ref. 9), (c. ) = 0.83 for 

l,max exp 

R n = 84 000 and (c, ) =1.55 for R n = 3 260 000. Linear interpolations 

c l ,max exp c 

between these values give (c, ) =0.84 for R n = 110 000 and (c, ) = 

l, max exp c ^ f max exp 

1.49 for R_ = 3 000 000. Therefore, (c, ) * 1.01(0.84/1.49) = 0.57 for 

c l, max exp 

NACA 642-015 with roughness for R c = 110 000. Since the lift-curve slope is not 

sensitive to Reynolds number, it is estimated to be (c. ) 0.110 (same as for 

1 a exp 

the higher Reynolds number). Therefore, F ~ 1.17 and a es ** 0.57/0.110 = 5°. The 
rolling-moment coefficient for this combination of parameters is below that for the 
higher Reynolds number and is more in line with the experimental result of refer- 
ence 2 (fig. 6) . 


Case B 

Velocity distributions of single and interacting trailing vortices and rolling 
moments on a following model measured in a wind tunnel are reported in reference 3. 
A measured velocity distribution for a merged vortex formed from two corotating 
equal-strength vortices (fig. 7) is used for case B since a corresponding measured 
C^ is given more accurately than others in reference 3. The velocity distribution 
is measured at x/B = 8.49 downstream of two semispan wings with B = 28.50 in. 

The rectangular following-wing model with an aspect ratio of 5.2 has an NACA 0012 
airfoil section, b = 4.56 in., c = 0.87 in., and R c * 50 000 in the experiment. 

For the NACA 0012 airfoil at the highest test Reynolds number (R c = 3 180 000) 

in reference 9, (c, ) = 0.103 and (c, ) = 1.56. Therefore, F = 0.103/ 

la exp l ,max exp 

0.094 = 1.10 and a es = 1.56/0.103 = 15° (a es = -15° for negative a). As shown 
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in figure 8, F increases from 0.130 to 0.143, but stall has no effect since 

the local incidence (arctan (V/UcJ ) does not exceed a es . (See fig. 7.) 

For R c = 43 000 (the lowest test Reynolds number in reference 9 for the 

NACA 0009 and NACA 0015 airfoils), (c. ) = 0.100 and (c n ) = 0.80 are 

la exp i r max exp 

average values from NACA 0009 and NACA 0015 data. Therefore, F = 0.100/0.094 = 1.06 
and a es = 0.80/0.100 = 8°. As shown in figure 8, F increases C^ from 0.130 to 

0.138, whereas a es = 8° significantly decreases it from 0.138 to 0.100, as might 
be expected from figure 7. 

Case B can be used to point out a potential problem with the interpretation of 
wind-tunnel results for low Reynolds numbers. The predicted C^ = 0.143 for 
R c = 3 180 000 (fig. 8) is much higher than the experimental result. On the other 
hand, the predicted C^ of 0.100 for R c = 43 000 (which is near the experimental 
R c = 50 000) is below the experimental result. These results indicate the rolling- 
moment coefficient for R c = 3 180 000 is about 40 percent higher than that for 
R c = 43 000. In other words, the rolling moment in flight may be much higher than 
that measured at a low Reynolds number in a wind tunnel. 


CONCLUDING REMARKS 

A method of computing rolling-moment coefficients for axial penetrations into 
a wake vortex is described in terms of simple modifications to a computer program 
which calculates wing-body aerodynamics. Wing twist is used to make allowances for 
the nonuniformity of the incoming airstream without considering the effects of the 
following wing on the vortex velocity distribution. Any computer program for 
subsonic wing aerodynamics for uniform free-stream flow can be modified in a 
similar way, provided that it has the capability of treating an arbitrary twist 
distribution across the full span of the wing. A program which treats arbitrary 
twist on an input semispan wing and reflects it across a plane of symmetry can also 
be used, provided that the input wing geometry can be separated arbitrarily far 
away from any reflected geometry. A section lift factor and an effective stall 
angle are used to make corrections for viscosity as well as for airfoil camber 
and thickness. 

The results of this method can be useful in studying the relative hazards 
of wakes from various vortex-generating aircraft configurations and in assessing 
and understanding the effect of stall on wind-tunnel results. Measurements of 
rolling moment for low Reynolds numbers can be misleading if stall occurs on the 
following wing. Flight data for higher Reynolds numbers will be less influenced 
by flow separation and may result in higher rolling moments. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
August 1, 1985 
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(a) Strong aerodynamic interference with no separation. 


y, in. 



(b) Weak aerodynamic interference with separation of 1 span. 


y, in. 



(c) Negligible aerodynamic interference with separation of 10 spans. 

Figure 2.- Effect of separation distance on interference between input 

and reflected wings (top view). 
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(a) Original numbering system for wing panels* 
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(b) Panels renumbered into left and right groups. 
Figure 3.- Numbering system of panels on input wing. 
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Figure 4.- Section lift curve measured in wind tunnel for NACA 0012 airfoil 

at R„ = 170 000 (ref. 9). 
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